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ABSTRACT 

Aims. We investigate the interaction of magneto- acoustic waves with magnetic network elements with the aim of finding possible 
signatures of the magnetic shadow phenomenon in the vicinity of network elements. 

Methods. We carried out three-dimensional numerical simulations of magneto- acoustic wave propagation in a model solar atmosphere 
that is threaded by a complexly structured magnetic field, resembling that of a typical magnetic network element and of internetwork 
regions. High-frequency waves of 10 mHz are excited at the bottom of the simulation domain. On their way through the upper 
convection zone and through the photosphere and the chromosphere they become perturbed, refracted, and converted into diff'erent 
mode types. We applied a standard Fourier analysis to produce oscillatory power-maps of the line-of- sight velocity. 
Results. In the power maps of the upper photosphere and the lower chromosphere, we clearly see the magnetic shadow: a seam of 
suppressed power surrounding the magnetic network elements. We demonstrate that this shadow is linked to the mode conversion 
process and that power maps at these height levels show the signature of three diff'erent magneto- acoustic wave modes. 

Key words. Sun: chromosphere - Sun: magnetic magnetic topology - Sun: oscillations - Magnetic Fields - Magnetohydrodynamics 
(MHD) 



1. Introduction 

The feature-rich solar chromosphere is one of the most disputed 
topics in present solar research (Rutten 2012). Debates are fu- 
eled by the ever-improving quality of the observations. For in- 
stance, high spatial resolution and high temporal cadence obser- 
vations reveal that oscillations in intensity and Doppler veloc- 
ities of chromospheric lines are suppressed in the vicinity, but 
not, or less, directly above quiet-Sun network elements. Judge 
|et al. ( |2001) , using SUMER UV spectral line time- series and 
TRACE continuum observations, referred to this phenomenon as 
magnetic shadow "since the network element appears to 'cast a 
shadow' over the neighboring ... internetwork region". These au- 
thors offered several explanations for the formation of the mag- 
netic shadow, favoring mode conversion of MHD waves if the 
plasma-yS (ratio of thermal to magnetic pressure) was about unity 
at the formation height of the shadow. 

In a subsequent study, [Mcintosh & Judge] ( |2001| ) used po- 
tential field extrapolations to reconstruct the magnetic field of 
the observed solar region, from which they concluded that the 
magnetic shadow is caused by closed magnetic field arches and 
not by mode conversion of MHD waves, and that the 13=1 
layer is likely located at the heights from which the UV con- 
tinuum emerges. [Krijger et"ar| ( |200l| also argued that at the ob- 
served height, mode conversion is unlikely to cause the magnetic 
shadow because this process would only become important at 
greater heights. 

Observations by |Vecchio et al. (20071 of the chromo- 
spheric spectral line of Can 854.2 nm, acquired with the 
Interferometric Bidimensional Spectrometer (IBIS) ( [Cavallini 
2006), most clearly revealed the ring and shadow-like appear- 
ance of the suppressed oscillatory power around quiet-Sun mag- 
netic network patches and at frequencies above the acoustic cut- 



oflT frequency. More recently, [Kontogiannis et al.| ( [2010a| ) showed 
similar results from Ha profiles obtained with the Dutch Open 
Telescope, revealing 'the discrete role of the magnetic field ... 
which guides or suppresses the oscillations, depending on its in- 
clination'. Kontogiannis et al. ( 2010b ) complemented this study 
with a potential magnetic field extrapolation and concluded that 
the magnetic shadow is associated with the < 1 environment. 

To shed light on the formation mechanism of the magnetic 
shadow, we applied detailed three-dimensional MHD simula- 
tions of magneto-acoustic waves to a model that shows both 
magnetic network and internetwork regions, although on a con- 
siderably smaller scale than the observations. We followed the 
simulation approach of Nu tto et al.| ([2012), which allows for 
the visualiziation of magneto- acoustic waves as they propa- 
gate across a complexly structured, dynamic solar atmospheric 
model. For the terminology used in this paper, we refer the 
reader to the introductory section of [Nutto et al. ( 2012| ). 



2. The background model 

The three-dimensional simulations were carried out with the 
CO^BOLD-code (Freytag et al.,2012 ), which solves the magne- 
tohydrodynamic equations for a fully compressible gas includ- 
ing radiative transfer and a realistic equation of state. 



2.1. Initial magneto-atmosphere 

The employed three-dimensional model extends from the sur- 
face layers of the convection zone throughout the photosphere 
to the middle layers of the magnetically structured chromo- 
sphere, corresponding to a height range from z = -1420 km 
to z = 1360 km, where z = km is defined by the height of 
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Fig. 1. Initial state. Left: emerging radiative intensity at the top 
boundary. Right: magnetic field strength on the r = 1 surface 
of the initial model. The horizontal line indicates the position of 
the vertical slice shown in Fig. 3] The contours in both panels 
display the equipartition level where Cs = ca- 
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mean optical depth unity, (r) = 1 . Rosseland mean opacities are 
used for the radiation transfer. The horizontal outline of the box 
is quadratic with the dimensions 4800 km by 4800 km. 

The numerical grid consists of 120 cells in each horizontal 
direction. In the vertical direction, the grid-cell size gradually 
decreases from 50 km near the bottom to 20 km in the atmo- 
spheric layers. Periodic boundary conditions apply to each hori- 
zontal direction, while the bottom and upper boundary are trans- 
mittent, allowing acoustic waves to leave the box without major 
reflection. The boundary conditions force the magnetic field to 
be vertical at the top and bottom boundaries but allow it to move 
freely in the horizontal direction. 

The left panel of Fig. [T] shows the bolometric intensity of 
the radiation that leaves the box through the top boundary in 
the vertical direction. It shows the regular granulation pattern 
of the solar surface and the intergranular lanes. The right panel 
displays the absolute magnetic field strength of the initial model 
at the surface layer of r = 1 . 

This initial background model was derived from a simula- 
tion that started with a vertically oriented, homogenous magnetic 
field of a strength of \B\ = 200 G. Because of the advective mo- 
tion of the granular cells, the magnetic field is concentrated in 
the intergranular lanes. In the course of the simulation, a strong 
flux concentration accumulates in the upper right quadrant of the 
box, at (x,y) = (3800, 3100) km. Here, field strengths of up to 
l^maxl = 1600 G are reached inside the flux tube at the surface 
of T = 1 . The left half of the box contains a rather weak, hori- 
zontally oriented magnetic field that resembles an internetwork 
region. 

2.2. Properties of the time-averaged model 

Starting from the initial model, the simulations were advanced 
for 1250 s. Figure [2] shows the time-averaged (1250 s) magnetic 
field for diff'erent optical depth levels, which is needed for the 
interpretation of the power maps that derive from the full length 
of the time series. 

The top panel of Fig. [2] displays the time-averaged absolute 
magnetic field strength, on the surface of the time-averaged 
position of optical depth r = 1 . Compared to the map of the ini- 
tial model shown in Fig.[T] the time-averaged magnetic field, es- 
pecially the strong flux concentration, appears difl'use because of 
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Fig. 2. Time-averaged absolute magnetic field strength, \Bl at 
the optical depth levels (a) ti = 1, (b) T2 = 0.0008, and (c) 
T3 = 6.7 X 10"^. The average is taken over the entire time se- 
ries of 1250 s. The horizontally averaged position of the consid- 
ered optical depth levels corresponds to the geometrical heights, 
(zin)} = km, (z(t2)) ~ 400 km, and (^(ts)) ~ 600 km. 
The white contour represents the time-averaged position of the 
equipartition level where ca = Cg. 

its continuous advective motion. Nonetheless, the time average 
still shows several strong magnetic field concentrations of about 
1600 G. At the transition to chromospheric levels, at T2 and T3, 
the magnetic canopies of several magnetic field concentrations 
merge and form a magnetically dominated chromosphere (cf. 
Sola nki & Steiner 1990 ). Below the magnetic canopy, the time- 
averaged background model contains a fairly weak, horizontally 
oriented magnetic field, similar to the initial model. 



3. Wave propagation 

Below, we show snapshots of the wave propagation simulations. 
We use the subtraction method described by Steiner et al. ( 2 007| ) 
for the visualization in which physical quantities of two simula- 
tion runs with and without an artificially introduced perturbation 
are subt racted from each ot her. Test results for this approach are 
given in [Nut to etaT]p010| ). 

The wave excitation is analogous to the two-dimensional 
case presented in Nutto et al. ( 2012 ), where plane-parallel waves 
are excited at the lower boundary. The same was performed for 
the three-dimensional case, except that the wave was now ex- 
cited along the entire plane of the lower boundary. The applied 
driver excites monochromatic plane-parallel waves with a fre- 
quency of 10 mHz. The waves are continuously excited for the 
whole duration of the simulation of 1250 s. The simulation can 
only be advanced for a limited amount of time after which the 
waves become swamped in the background of the residual large- 
scale velocities (see [Nutto et al.|20f2 ). Therefore, we choose the 
relatively high frequency of /o = 10 mHz to increase the reliabil- 
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Fig. 3. Time instances of the wave-propagation simulations 
showing the logarithm of the absolute velocity perturbation, ^|v|. 
Time increases from top to bottom. The vertical slice is taken 
along the xz-plane at the position of a strong magnetic flux tube 
(see horizontal line in the right panel of Fig.[T]). The green con- 
tour corresponds to optical depth r = I. The white contour dis- 
plays the equipartition level. 



ity of the analysis for the available limited time period and also 
for convenient visualization of the waves. One can expect a sim- 
ilar behavior of propagating waves for all frequencies above the 
acoustic cut-ofl' frequency ( Khomenko & Collados 2006), there- 
fore waves with /q = 10 mHz should be a good representation 
for all waves above the acoustic cut-off' frequency of 5 mHz. 

At the bottom of the simulation domain and throughout the 
simulated portion of the convection zone, we have » 1 so that 
the excited wave is a fast predominantly acoustic mode. The up- 
per panel of Fig.[3]shows the wave as it reaches the solar surface. 
Mode conversion has not significantly taken place yet, because 
in the center of the flux tube, the angle between the wave vector 
and the magnetic field (attack angle) favors the transmission of 
the acoustic mode ( [Cally|20Q7 ). 

In the lower panel of Fig. |3| the wave has traveled through 
the photosphere and has reached the lower layers of the chro- 
mosphere. Clearly, the wave field shows a diff'erent propaga- 
tion behavior inside and outside the magnetic canopy of the flux 
concentration. Below the magnetic canopy, and hence between 
the magnetic network elements, the wave propagation is nearly 
undisturbed. For regions inside the magnetic network elements, 
viz., above the equipartition leveQ the wave fronts are vertically 
directed and belong to the fast magnetic mode that is excited be- 
cause of mode conversion at the equipartition level. Owing to 
steep gradients of the Alfven speed in the vertical and horizon- 
tal directions, the fast magnetic mode is strongly refracted. In 
the apex of the refractive path, the fast magnetic mode travels 
horizontally with the wave front being aligned with the vertical 
direction. 

In the very center of the magnetic flux element and coexis- 
tent with the vertical wave fronts of the fast magnetic mode, we 



^ The equipartition level is defined as the location where the Alfven 
and the sound speed are of equal magnitude, ca = Cs 



see the transmitted slow acoustic mode. Its wave front keeps the 
horizontal orientation from below and travels vertically upward 
along the magnetic field lines of the flux concentration, reaching 
the higher layers of the solar atmosphere much later than the fast 
mode. 



4. Power maps 

We restricted our investigation to power maps of the vertical ve- 
locity perturbations, Sv^, because these can be directly compared 
to power maps of line-of-sight velocities observed at the center 
of the solar disk. 

The power maps were taken at constant optical depth levels 
n = 1, T2 = 0.0008, and T3 = 6.7 x 10"^ in the atmosphere. The 
time-averaged positions of these optical depth levels correspond 
to the geometrical heights (z(ti)) = km, (z(t2)) ~ 400 km, and 
<z(t3)) ^ 600 km. The time series was sampled with a cadence of 
At = Is. The power of observed Doppler shifts of active regions 
are usually normalized to regions of the quiet Sun that show the 
least magnetic activity. In our simulations, where the whole sim- 
ulation domain is threaded by a magnetic field, we cannot deter- 
mine a corresponding reference value. Therefore, we normalized 
the power to the value of the spatially averaged power of each 
optical depth level. 

The derived power maps for the upper photosphere, T2, and 
the lower chromosphere, T3, presented in Fig. |4j yield rings of 
suppressed power (indicated by ellipses) around strong magnetic 
field elements. These resemble very much the magnetic shadows 
observed by Vec chio et al.| p007j ) in the lower chromosphere 
above magnetic network elements. Three separate regions of dif- 
ferent oscillatory strength can be identified. The strongest oscil- 
latory power occurs in the internetwork regions (areas inside the 
equipartition contour where jS » 1), while intermediate power 
can be found in the very center of the magnetic network elements 
(inner ellipses). The suppressed oscillatory power, and hence the 
shadow-like appearance, occurs in the regions between the inner 
and outer ellipses. 

The variable strength of the oscillatory power for the diff'er- 
ent regions is caused by the presence of various modes of the 
magneto-acoustic wave. In the internetwork regions the strong 
power results from the nearly undisturbed fast acoustic mode 
which propagates upward through the almost field-free medium 
(cf. Fig. |3]). Here, the wave mode primarily possesses a longi- 
tudinal polarization with velocity perturbations, Sv^, along the 
propagation path and line-of-sight. Thus, the full strength of the 
oscillation shows up in the power map. 

In the very center of the magnetic flux tubes, where the wave 
vector and the magnetic field vector are parallel, the fast acous- 
tic mode has been primarily transmitted to the slow acoustic 
mode. This upwardly propagating slow acoustic mode of lon- 
gitudinal polarization causes the intermediate oscillatory power 
in the very center of the magnetic network elements. Compared 
to the power of the acoustic mode in the internetwork region, 
the reduction of power is owed to the mode conversion process, 
where a fraction of the acoustic energy has been converted into 
the fast magnetic mode. 

As shown in Sect. [3] the fast magnetic mode propagates in 
the region between the internetwork and the center of the net- 
work elements. The suppression of the oscillatory power can 
now be understood in terms of wave polarization — specifically 
by the orientation of the wave vector with respect to the mag- 
netic field and line-of-sight. In the apex of the refractive wave 
path, the fast magnetic mode travels basically perpendicular to 
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Fig. 4. Power maps of the vertical velocity perturbations, Sv^, taken at the photospheric optical depth levels T3 = 0.0008 (a), and at 
the lower chromospheric layer, T4 = 6.7 x 10"^ (b). The white contours show the time- averaged position of the equipartition level. 
The ellipses mark regions where the magnetic shado w can be identified . Between the large and the small ellipses, the characteristic 
ring of suppressed power is apparent (cf. Fig. 2 of ,Vecchio et al.| ( [2QQ7| )). 



the magnetic field with longitudinal wave polarization and con- 
sequently, with velocity perturbations in the direction of the hor- 
izontally directed wave vector. Therefore, along the propagation 
path, the initial vertical velocity perturbations, Sv^, are gradu- 
ally transformed into horizontal velocity perturbations, Svx,y, and 
back again. Thus, the oscillatory power of the vertical velocity 
perturbations, Sv^, corresponding to line-of-sight velocities at the 
solar disk center, are suppressed around magnetic network ele- 
ments, causing th e appearance of a magnetic sh adow. 

According to |Khomenko & Cally (2012 and references 
therein), part of the fast magnetic mode can convert into an 
Alfven wave, which then propagates along the magnetic field 
lines. The fast- to- Alfven conversion takes place most efficiently 
in the evanescent tail of the fast wave beyond the apex of the re- 
fractive wave path (reflection point), w hich extends t hroughout 
the upper chromosphere and beyond ( Kho menko & C ally 2012 ). 
Therefore it cannot be expected to leave imprints on the power 
maps of Fig.|4j which pertain to the low chromosphere and upper 
photosphere, although it may become detectable in the core of 
Ha. Moreover, in our case, the apex of the refractive wave path 
occurs in regions threaded by mainly vertically directed mag- 
netic fields. Hence, the excited Alfven waves would not signifi- 
cantly contribute to line-of-sight velocity perturbations because 
of their transversal character. Consequently, we can consider the 
generation and propagation of the fast magnetic mode as the pri- 
mary cause for the appearance of the magnetic shadow. 



5. Conclusion 

Power maps of high-frequency waves observed in chromo- 
spheric spectral lines show a reduction in power in the surround- 
ings of network magnetic patches reminiscent of a shadow. The 
responsible physical process for this appearance has remained 
ambiguous so far. With three-dimensional radiation magneto- 
hydrodynamic simulation, we can now reproduce the magnetic 
shadow phenomenon and we conclude that it is tied to the pro- 



cess of mode conversion. The excitation of the fast magnetic 
mode in the expanding region of the network magnetic field and 
its refractive wave path reduces the oscillatory power of the ver- 
tical velocity perturbations, Sv^, in the surroundings of magnetic 
network elements. 

Furthermore, we can show that the power maps taken in the 
upper photosphere and lower chromosphere show the oscillatory 
signatures of diff'erent magneto-acoustic wave modes in three 
diff'erent phenomenological regions: the fast acoustic mode in 
the internetwork, the fast magnetic mode with possible superpo- 
sition of Alfven waves in the magnetic canopy region, and the 
slow acoustic mode in the center of magnetic network elements. 
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